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The doping of quantum Heisenberg antiferromagnets on the kagome lattice by non-magnetic impurities is 
investigated within the framework of a generalized quantum dimer model (QDM) describing a) the valence 
bond crystal (VBC), b) the dimer liquid and c) the critical region on equal footing. Following the approach 
by Ralko et al. [Phys. Rev. Lett. 101, 117204 (2008)] for the square and triangular lattices, we introduce 
the (minimal) extension of the QDM on the Kagome lattice to account for spontaneous creation of mobile 
S=l/2 spinons at finite magnetic field. Modulations of the dimer density (at zero or finite magnetic field) 
and of the local field-induced magnetization in the vicinity of impurities are computed using Lanczos Exact 
Diagonalization techniques on small clusters (48 and 75 sites). The VBC is clearly revealed from its pinning by 
impurities, while, in the dimer liquid, crystallization around impurities involves only two neighboring dimers. 
We also find that a next-nearest-neighbor /erromagneftc coupling favors VBC order. Unexpectedly, a small size 
spinon-impurity bound state appears in some region of the the dimer liquid phase. In contrast, in the VBC phase 
the spinon delocalizes within a large region around the impurity, revealing the weakness of the VBC confining 
potential. Lastly, we observe that an impurity concentration as small as 4% enhances dimerization substantially. 
These results are confronted to the Valence Bond Glass scenario [R.R.R Singh, Phys. Rev. Lett. 104, 177203 
(2010)] and implications to the interpretation of the Nuclear Magnetic Resonance spectra of the Herbertsmithite 
compound are outlined. 



I. INTRODUCTION 

The spin- 1/2 quantum Heisenberg antiferromagnet (QHAF) 
on the kagome lattice (see Fig. [TJa)) is believed to be the 
paradigm of frustrated quantum magnetism. It is also the 
best candidate to observe exotic (non-magnetically ordered) 
phases predicted by theorists. Among those, an algebraic spin 
liquid', a gapped dimer -or spin- liquicPS or a number of 
valence bond crystals (VBC) with 12-sit^, 18-sitd^and 36- 
site^ ^ unit cells or with columnar dimer ordeJ^have all been 
proposed as possible ground state (GS) of the kagome QHAF. 

On the experimental side, Herbertsmithite^ is an almost per- 
fect realization of a kagome QHAF. In this material Copper 
atoms carrying spin- 1/2 degrees of freedom (oxidation state 
+2) are located on the sites of (very weakly coupled) kagome 
layers. In addition to the crystallographic structure, the close- 
ness to an ideal system is supported by the remarkable ab- 
sence of any magnetic ordering at the lowest attainable tem- 
peratureJ^, as expected for a highly frustrated quantum mag- 
net. 

However, small deviations from a perfect kagome QHAF 
are known in this material. In addition to a small 
Dzyaloshinskii-Moriya (DM) anisotropy", Herbertsmithite 
contains a small amount of non-magnetic Zinc impurities sub- 
stituting Copper atoms in the kagome planes. Nuclear Mag- 
netic Resonance (NMR) offers a fantastic tool to probe im- 
purities via the change induced by a nearby impurity on the 
local magnetic field at the nucleus whose resonance is being 
measured. The first estimations from NMR'^ of the Zinc con- 
centration range from 6 to 10 %. It is not clear therefore, 
which experimental features in Herbertsmithite can really be 
attributed to intrinsic properties of the kagome QHAF and 
which of them are induced e.g. by Zinc substitution. 

One of the strong motivations of this work is to investi- 
gate the effect of spinless impurities in quantum disordered 




FIG. 1: (Color online), (a) Kagome lattice consisting of a hexago- 
nal lattice of comer-sharing frustrated triangular units. The Heisen- 
berg AFM interaction acts on all bonds of every triangle, (b) A non- 
magnetic static impurity (lozenge) in the kagome lattice is described 
as a vacancy. Mobile spinons (arrows) carry spin- 1/2 degrees of free- 
dom. (*) The N — 75 cluster contains only a single impurity and no 
spinon. (c) The Mahattan distance d is obtained by counting the num- 
ber of unit segments (oriented along the two crystallographic axis) on 
the path joining the impurity to a given bond center. Non-equivalent 
bonds are labelled by greek letters, (d) A spinon (carrying spin- 1/2 
shown by a blue arrow) on a triangle cannot hop on the sites of the 
same triangle, (e-f) A spinon can hop at further distances leading to a 
120-degree (c) or 180-degree (back-flow) rotation (d) of the involved 
dimer. 



(i.e. non-magnetic) phases which are the most serious candi- 
dates for the Herbersmithite ground state (GS). Starting from 
a simple description of this material in terms of a SU(2)- 
symmetric kagome spin- 1/2 QHAF (neglecting the small DM 
anisotropy) recent calculations of NMR spectra'^ seem to re- 
flect key experimental features'^. On the other hand, the lo- 
cal susceptibilities (Knight shifts) obtained by high tempera- 
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ture series expansiorP^ were claimed to agree with an alge- 
braic (gapless) spin liquid phase. Here we provide an alter- 
native analysis based on a generalized quantum dimer model 
(QDM) which can naturally describe a) topological (so-called 
Z2) spin-liquid, b) VBC and c) the critical region on equal 
footing and distinguish features around impurities which are 
specific to each phase. The paper is organized as follow: first, 
in Section|Ilj we describe the QDM framework, the modeling 
of the static dopants and the simplest extension of the QDM 
to account for a small magnetization (induced by a finite mag- 
netic field) in the ground state. In Section III the results for 
a single imbedded impurity are described. This includes the 
investigation of the pinning of the VBC and the calculation of 
the spin density of a doped mobile spinon around the vacancy. 



Then, in Section IV we move on to the case of several impuri- 
ties and show how the dimerization pattern can be enhanced. 
A discussion of the experimental NMR data in the light of 
our results is given in Section [V] as well as some concluding 
remarks in SectionlVTl 



II. THEORETICAL DESCRIPTION 



to -1/4) has also been introducecP! Following Ref. |22l we 
shall consider a linear interpolation between the two Hamil- 
tonians, i/^ff (A) — Ai?Hois + (1 ^ '^)^rk- As we shall see 
next, it is remarkable that this "generalized" QDM can de- 
scribe topological dimer-liquid, VBC and the critical region 
on equal footing. Whether this model (for A = 1 or around 
this value) describes faithfully the original (unprojected) mi- 
croscopic QHAF on the Kagome lattice is still under debate. 
However, lots of insights on how impurities behave generi- 
cally in the quantum disordered phases described by this ef- 
fective model can be obtained. 



^eff = ^6 0+ (0+Z^+Z7) 

+^10 + ^12!^ 



A. Impurities described as vacancies 

Let us first discuss qualitatively the effect of a single im- 
purit}G4l. 2inc impurity is a spinless ion and therefore can 
be described as a vacant site as depicted in Figs. [T|b,c). In 
other words, the bond exchange interaction Si • (in units 
of the exchange coupling J) acts on all bonds of the kagome 
lattice except on the four bonds connected to the impurity. Us- 
ing Lanczos exact diagonalizations (LED) of small clusters of 
the kagome QHAF it has been shown that a single impurity 
tends to 'localize' two singlet bonds next to iiMlSl Although 
this local phenomenon is well captured even on very small 
clusters'^, LED of the QHAF do not allow to investigate reli- 
ably the perturbation of the media at larger distances from the 
impurity, which can be probed e.g. by NMR. Therefore, we 
shall use here an effective description based on the re cently 
developed QDM allowing to handle larger clusterJSEo] 



B. Competing singlet phases in the pure system 

Starting from a projection of the QHAF on the (non- 
orthogonal) NN singlet basis^' (an approximate resonating 
VB description) followed by a transformation onto an ad-hoc 
orthogonal 'dimer' basis and, finally, by a series expansion 
(in terms of a small overlap parameter a = l/-\/2), one gets 
an effective QDM of the type shown in Fig. |2] Here we use 
the ampUtudes obtained in Ref.|20l Vq = 1/5, Vs = 2/63, 
Vio = 1/255, V12 = 0, Je = -4/5, Jg = 16/63, Jw = 
— 16/255, J12 = (in units of J) which defines our "Heisen- 
berg" effective QDM i/ncis- For convenience, the labels 7 
in the amplitudes and Vy correspond to the lengths of the 
associated resonance loops. An exactly solvable "Rokhsar- 
Kivelson" (RK) QDM i?RK defined by vanishing diagonal 
amplitudes — and equal kinetic amplitudes (set here 



FIG. 2: The effective QDM used here. A sum over all the hexagons 
of the lattice is implicit. Kinetic terms (J^) promote cyclic permu- 
tations of the dimers around the loops (shown as thick lines) and 
diagonal terms (V^) count the numbers of "flippable" loops. 

So, prior to the introduction of impurities, let us give a 
brief summary of the results of the pure generalized QDM 
(in zero magnetic field). Generically, short or long-range or- 
deiing of the dimers can happen yielding to dimer liquids or 
VBC, respectively. Numerical results^ show that a 36-site 
unit cell VBC (shown schematically in Fig.[3|a)) is stabilized 
for A > Ac — 0.94, in particular at the "Heisenberg" point 
A = 1. A remarkable (topological) gapped Z2 dimer liquid 
which bears an exact analytical form at the RK point A = is 
stable for A < Ac as shown in the phase diagram of Fig.[3jb). 
Besides, at exactly A = 1, J12 vanishes so that the undeter- 
mined chiralities of the pinwheel (or star) resonances of the 
VBC pattern (see e.g. Ref. '2T for details) lead to an extra 
(Ising-like) macroscopic degeneracy. These findings at A = 1 
are in perfect agreement with recent series expansions on the 
QHAF^. Note that the proximity to a Quantum Critical Point 
(QCP) suggests that dimer fluctuations remain strong and that 
the VBC amplitude is weak at A = 1. However it should 
be stressed here that the derivation of the A = 1 effective 
QDM is only semi-quantitative since it relies i) on the projec- 
tion of the Heisenberg Hamiltonian onto the singlet subspace 
spanned by the NN VB and ii) it is subject to a (controlled) 
truncation of a series expansion. Besides, the real material 
might contains longer range exchange interactions. This sug- 
gests that a faithful description of the material in terms of a 
generalized QDM might not correspond exactly to A = 1. We 
can however consider A as a "phenomenological" parameter 
which enables to "tune" the system in one of the two quan- 
tum disordered phases or in the critical region in between to 
investigate the generic role of impurities in each case. 
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FIG. 3: (a) Schematic representation of the 36-site VBC supercell 
(from Refs. I20l22t . (b) Phase diagram of the pure QDM. In the 
present study, one moves along the (blue) A-axis. 



C. Extending the QDM to finite magnetic field 

Spinon hoppings - In an NMR setup a small magnetic field 
is applied to slightly polarize the system. The local site- 
dependent magnetization is then experimentally accessible via 
the measured Knight shift. However, ground states with fi- 
nite magnetization cannot be addressed theoretically within 
the above QDM which only describes the singlet subspace dy- 
namics. Nevertheless, a finite field/magnetization setup can be 
realized theoretically by introducing extra spin- 1/2 degrees of 
freedom (named "spinons" hereafter) in the QDM description 
as first proposed in Ref. 27 Physically, such a spinon (polar- 
ized along the magnetic field) can be viewed as resulting from 
the breaking of a singlet bond by the introduction of a vacant 
site (impuiity). Spinon nearest-neighbor pairs can also ap- 
pear or disappear by exchanging with dimers (see below). The 
bond exchange interaction of the QHAF leads to the motion 
of the spinons followed by some dimer "backflow''^. Since, 
the generalized QDM has been derived in Ref. '20' assum- 
ing a fermion representation of the (original) SU(2) dimers, 
consistency implies that the mutual (bare) statistics of the 
spinons should be fermionic^'^. Here, we extend the proce- 
dure of Ref |20] to derive from the microscopic QHAF the 
simplest effective Hamiltonian (in lowest order) govern- 
ing the dynamics of the (fermionic) spinons sufficient to cap- 
ture the main features of spinon delocalization. The relevant 
processes are shown in Figs. [TJd-f). The overlap matrix ele- 
ments between the initial and final SU(2) configurations |(/)) 
and IV') involved in these processes are O0,V' = ^1/2- The 
corresponding Heisenberg matrix elements Ti,^,^ equal +1/2 
(in units of J) for the last two processes of Figs. [Tj6,f) and 
vanishes for the first one of Fig. [ijd). Note that, for consis- 
tency with the definition of the QDM of Fig. [2] we subtract a 
singlet-bond energy and include a 4/3 multiplicative factoi^^. 
Using the expansion scheme given by (the first two terms of) 
Eq. (41) in Ref. '20', we deduce easily that, in the effective 
QDM, the first process vanishes and the two remaining ki- 
netic processes shown in Figs.[TJe,f) involve the same "spinon 
hopping" amplitude t^-p ~ +1/2. This procedure also gen- 
erates (at order a^) a potential term of amplitude V^sp = 1/4 
which "counts" the number of possible "spinon hopping" of 
each kind. Stricly speaking, this (minimal) effective spinon 



Hamiltonian is most relevant to the physical case A = 1 but it 
is also interesting to investigate its properties in an extended 
range e.g. 0.5 < A < 1. 



(a) (h) 




h = hj^O 



FIG. 4: Periodic clusters used in this study; (a) The 5\/3 x 5\/3 
(A*' — 75) cluster with a single static impurity, (b) The 4\/3 x 4\/3 
(A'^ — 48) cluster doped with — 2 static impurities and A'sp = 
2 mobile spinons (polarized along the field). 

Spinon chemical potential - The density of spinons Ugp — 
Nsp/N is tuned by the magnetic field H. Indeed, A'sp spinons 
polarized along the field gain Zeeman energy —{h/2)Nsp, 
where h = ^giigH, which is balanced by a creation energy 
term (Ab/2) A'sp. Note that we rescale here the magnetic field 
by the same 4/3 factor as the QDM ^"^^ . Physically, the param- 
eter Ab corresponds to the energy cost of breaking a sin glet 
bond into a triplet made of two polarized static NN spinons^Sl 
Practically, {h — Ab)/2 plays the role of a chemical potential 
for the spinons and the (grand canonical) GS energy becomes, 

Eg^{h,Nin,p) = Mm/N^^{E^s[N^p,Nin,p]-^{h-A-B)N^p} , 

(1) 

where Ech[Nsp, Nimp] is the GS energy of the QDM depicted 
in Fig.|2]doped with iVsp mobile spinons and iVimp (random) 
impurities and where a minimization over A'sp has been car- 
ried out. Although, strictly speaking, one expects A^ = 1 
(in units of X/"^ within the procedure of Ref. 20, to derive 
the effective QDM, its "optimum" value for a given order of 
the expansion scheme is unclear. However, Ab can also for- 
mally be related to the physical spin gap Ag of the pure sys- 
tem, i.e. the minimal energy cost to create a triplet (S = 1) 
two-spinon state in zero magnetic field (see below). Lastly, 
we note that two spinons can form a (triplet) bound state^^ 
whenever Ecb[2, 0] + E^s[0, 0] - 2£;eff[l, 0] < (in the pure 
system). 



m. SINGLE IMPURITY RESULTS 

A. Local dimer modulation 

As a preliminary study, we now start by just inserting a sin- 
gle spinless impurity (no spinon). As stated above, we use the 
modelisation of the impurity in terms of a vacant site as shown 
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in Figs.[TJb,c). In the framework of the generalized QDM, the 
vacancy suppresses the resonant and diagonal terms of Fig.|2] 
whose loop contains the vacancy sitd^I 
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FIG. 5: (color online) (a) Spectrum versus A of the effective QDM 
on a N=75 site periodic cluster containing a single impurity and with 
zero magnetization. The 8 lowest energy levels are shown in each of 
the 4 symmetry sectors. Inset: generalized susceptibility (see text) 
showing the location of the phase transition. For A < Ac, the large 
spectral gap above the quasi-degenerate topological GS correspond 
to the singJe-vison gap of the Z2 liquid, (b) Average dimer density 
versus Manhattan distance from the impurity on the same N=75 clus- 
ter. The short distance behavior is separated in the inset for clarity. 
The mean value corresponds to 1/4 dimer per bond. 

In the very diluted impurity limit, one does not expect 
to fundamentally modify the bulk properties of the quantum 
magnet. However, a single impurity can e.g. pin VBC order 
and, in some way, offers a very interesting probe of its exis- 
tence. To test this, we have considered the 5\/3x 5\/3 periodic 
cluster of Fig. [4|a) with N ^ 75 sites and a single impurity 
(to comply with the hard-core dimer constraints such a cluster 
has to contain an odd number of sites j^S The full spectrum of 
this cluster is shown in Fig.|5ja) as a function of the parameter 
A. Here we have distinguished between the four topological 
sectors (for such a torus topology)^^ and used the reflection 
symmetry w.rt. one of the (72 -axis passing through the impu- 
rity site (cr = ±1). A quantum phase transition^^ occurs as a 
function of A for A = Ac — 0.9373 as evidence from the sharp 
peak in the GS generalized susceptibility x\ = d^EQs/dX^ 
shown in the inset^. A topological Z2 liquid (with four-fold 
degeneracy on a torus) is stable for A < Ac- At A = Ac the 
single-vortex (or "vison") gap"*^ corresponding to the first odd 
(cr = — 1) energy excitation vanishes. For A > Ac the system 
spontaneously breaks reflection symmetry (the GS is two-fold 
degenerate with cr = ±1 quantum numbers) as expected in the 
36-sites VBC. These results and, in particular, the estimation 
of the critical value Ac — 0.9373 are fully consistent with the 
results obtained on a pure 108-sites clusteil^S. 

It is interesting to examine the behavior of the dimer den- 
sity as a function of the Manhattan distance (see Fig. [ijc)) 
from the impurity, as reported in Fig. |5jb) for different val- 



ues of A between 0.5 and 1 . One striking feature is the large 
value of the dimer density on the two bonds facing the impu- 
rity (distance d = 2). This can be interpreted as a (singlet) 
dimer crystallization as found in prior studies of the Heisen- 
berg quantum antiferromagnel^— . Our study reveals that this 
phenomena is very robust and depends weakly on the sup- 
posed phase, liquid (A < Ac) or crystalline (A > Ac), of 
the model. Therefore this feature cannot be used practically 
as an experimental fingerprint. In contrast, when considering 
longer distances, strong differences in the bond modulation 
occur in the supposed liquid and VBC phases. For A < Ac, 
in the Z2 dimer liquid, the dimer density becomes very uni- 
form beyond distance 1 . When increasing A the appearance of 
strong modulations is exactly correlated to the crossing of the 
critical value and the entering in the VBC phase. For A = 1, 
the relative bond modulation amplitude is of order 20%. 



B. Effect of a next-nearest neighbor exchange coupling 

Recent calculations using variational Gutzwiller-projected 
wavefunctions have shown that a small next-nearest neighbor 
(NNN) ferromagnetic coupling j2 (measured in units of the 
NN exchange) may favor VBC ordering^^. A ferromagnetic 
coupling should then enhanced the dimer response around 
impurities, an effect which can be tested within the QDM 
framework by considering again the N = 75 cluster doped 
with a single impurity (no spinon). From Ref. 20, the new 
QDM amplitudes and V^, 7 = 6, 8, 10 have simply to 
be multiplied by extra coefficients (1 + A^) which depend 
on j2- is given by the ratio of the corresponding QHAF 
Hamiltonian matrix elements and, from Eq. (12) of Ref. |20l 
giving the general form of these matrix elements, one gets 

A = h E(.j)e-y,NNN e^j/(7/2 + E(.,j)e7,NN e^j) where each 
loop topology has now to be distinguished and j — — 1 (+1) 
if the sites i and j on the loop 7 are separated by an odd 
(even) distance. Note that for 7=12 (star resonance), the 
NN QHAF matrix element 7/2 + j)£-y nn ^iJ vanishes so 
that the leading contributions to J12 and V12 are directly pro- 
portional to j2- Using the general expansion scheme given by 
Eq. (41) of Ref. l20lit is straightforward to obtain the leading 
contribution (at order a^^ and a^" respectively), J12 = — §^2 

andl/12 --2|;J2. 

The excitation spectrum for A = 1 shown as a function 
of j2 in Fig. |6|a) reveals two gapped phases separated by a 
narrow gapless region located around ^ 0.05 (which might 
become a quantum critical point in the thermodynamic limit). 
The j2 < ^2 region (including j2 — 0) corresponds to the 
above pinned VBC as confirmed by the GS quasi-degeneracy 
and the enhanced VBC patterns shown in Fig. |6|b) for weak 
ferromagnetic j2- In contrast, for AFM NNN exchange j2 > 
^2, only a small dimerization is observed in Fig.|6|b), except at 
short distances from the impurity. This might be the signature 
of a dimer liquid or a very weakly pinned VBC and more work 
(on pure systems) is needed to resolve this case. 
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FIG. 6: (color online) (a) Spectrum versus j2 (NNN coupling in units 
of the NN coupling) of the A = 1 effective QDM on a N=75 site pe- 
riodic cluster containing a single impurity and with zero magnetiza- 
tion. The lowest energy levels are shown in each of the 4 symmetry 
sectors, (b) Average dimer density versus Manhattan distance from 
the impurity on the same N=75 cluster. The short distance behavior 
is separated in the inset for clarity. The mean value corresponds to 
1/4 dimer per bond. 



C. Does a spinon bind to an impurity ? 

Notoriously, spinons in dimer liquids (VBC) are 
known to be deconfined (confined). However, confine- 
ment/deconfinement is an asymptotic (e.g. long-distance) 
property of the system and unexpected behaviors can appear 
at shorter distance (comparable e.g. to the typical impurity 
separation). Indeed, a spinon could bind to an impurity 
(acting here as a "static spinon") even in a dimer liquid 
or, reversely, a spinon could spread over a very large area 
(e.g. beyond the cluster sizes available numerically) even 
in a VBC. To investigate such issues, we now consider a 
4a/3 X 4^3 (48-sites) cluster^"* (j2 = 0) doped with a 
static vacancy and a mobile spinon. The maps of both the 
(bond) dimer density and the site magnetization are shown 
in Fig. [T] for three distinct regimes : (i) A = 0.5 < Ac, 
deep in the dimer liquid phase; (ii) A = 0.9 < Ac and 
(iii) A — 1 > Ac, in the dimer liquid and VBC phases 
respectively, in the vicinity of Xq- The dimer density maps 
in this lightly polarized system closely resemble the previous 
zero field data of Fig.jSj^b) : the dimer liquid phase at A = 0.5 
exhibits a very homogeneous dimer density apart from the 
two "crystallized" bonds next to the impurity while strong 
dimer patterns are clearly visible in the VBC phase. Note that 
the dimer pattern gets more pronounced in the dimer liquid 
phase when approaching the dimer liquid- VBC transition, as 
e.g. for A = 0.9. The local magnetization maps of Figs.|7|b) 
& (c) are entirely correlated with the dimer density maps: 
for A = 0.9 (liquid) the magnetization is slightly depressed 
on the four bonds NN to the impurity and very uniformly 
distributed beyond. For A = 1 (VBC) a strong modulation 
in the local magnetization appears. In other words, the 




a) A=0.5 



b) A=0.9 



c) A=1.0 



FIG. 7: (color online) Real-space maps of the the dimer density 
(segments) and site magnetization (circles) computed on a periodic 
A'^ — 48-site cluster doped with a single impurity and a mobile 
spinon, for different values of A (as shown on the plots). The radii 
(width) of the circles (segments) are directly proportional to the val- 
ues of the corresponding observables. 



magnitudes of the dimer and of the (magnetic field induced) 
spin density modulations seem to scale with each other. The 
behavior in the dimer liquid phase is fully consistent with the 
expected spinon deconfinement mechanism found by LED 
of small clusters of the Kagome QHApfSEl In contrast, the 
VBC sample does not seem to show at the length scale of the 
cluster the (opposite) confinement behavior found e.g. in the 
VBC phase of the QHAF on the checkerboard lattice23. 

Although probably not directly relevant to experiments, the 
data of Fig. |7ja), deep in the dimer liquid phase, are of par- 
ticular theoretical interest : here, although the dimer pattern 
is very uniform (as expected in a dimer liquid phase) and no 
confining potential is expected, the spinon still forms a tight 
bound state with the impurity. We believe however that this 
behavior is not generic and depends strongly on the details of 
the spinon Hamiltonian : here the bound state is stabilized by 
an increase of spinon kinetic energy in the vicinity of the im- 
purity. Whether such bound state can appear spontaneously 
without magnetic field depends on the value of the spinon 
chemical potential i.e. on A^. If this is the case, the impurity- 
doped system becomes gapless for magnetic excitations (see 
below). 



IV. FINITE IMPURITY CONCENTRATION 



We now consider the case of A'imp — 2 impurities doped 
into the same cluster of = 48 sites, attempting to (crudely) 
mimic a finite density of impurities njmp — Nnnp/N of order 
4%. For simplicity, we fix the relative position between the 
two impurities at a "typical" distance, as shown in Fig. |4|b). 
We shall also consider the same values of the parameter A as 
above in the investigation of a single impurity, namely A = 
0.5, 0.9 and 1, for comparison. Hereafter, we also set j2 = 0. 



A. Enhanced dimerization 

Let us first assume zero total magnetization e.g. A^sp = 0. 
The corresponding data are shown in Figs.[8]^a-c). We clearly 
observe a sizable increase of the average dimerization com- 
pared to A^imp = 1- This is particularly clear for A — 0.9, 
a parameter corresponding, in the pure system, to the dimer 
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a) A=0.5 



b) A=0, 



c) A=1.0 



FIG. 8: (color online) Real-space maps of the dimer density (seg- 
ments) on a periodic TV = 48-site cluster doped with two impurities. 
Same conventions as Fig.|7] 
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liquid phase close to the phase transition at A = Ac. In fact, 
with two doped impurities, the new dimerization pattern in 
Fig.[8|b) resemble very closely the one at A = 1 in Fig.jSjc). 
It is interesting to also note that the local dimerization patterns 
around impurities can be different. Indeed, e.g. in Figs.|8jb) & 
(c), while one impurity is surrounded by two (quasi-isolated) 
strong dimer bonds, the second impurity stabilizes two neigh- 
boring symmetric resonating hexagons. All these findings are 
in fact consistent with the proposed picture of the Valence 
Bond Glass^^ 



FIG. 9: (color online) Spin gap of the pure (filled symbols) and 
two-impurity (empty symbols) doped 48-site cluster vs the param- 
eter Ab- 




a) A=0.5 



b) A=0.9 



c) A=1.0 



B. Spin gap 



FIG. 10: (color online) Real-space maps of the dimer density (seg- 
ments) and site magnetization (circles) computed on a periodic A'^ = 
48-site cluster doped with two impurities and two mobile spinons. 
Same conventions as Fig. [7] 



We now considerer the possibility that the ground state car- 
ries a small magnetization proportional to the (doped) spinon 



density m 



If a spin gap exists at /i = 0, we expect 



that a large enough magnetic field will eventually suppress the 
gap and induce a small magnetization. In our picture, the spin 
gap is the minimal energy to create a pair of mobile spinons. 
From Eq. ([T]i, the spin gap at low field h is therefore given by: 



As(ft.,nimp) 



Ecff [2, iVimp] 

As(0,niinp) - 



h 

(2) 



where Ab can be taken as a phenomenological parameter 
(considering the fact that our spinon hamiltonian is only a 
minimal model and does not include the potentially important 
higher-order longer-range hoppings). As in the conventional 
picture, the spin gap (if any) vanishes linearly at a critical 
field. We have numerically computed the spin gap in the same 
iV = 48 cluster doped with two impurities (same A's also as 



above) i.e. for an impurity concentration rii, 



Mmp/A^ 



4%, and compare it to the equivalent data in the pure system 
(nimp = 0) in Fig. [9] Interestingly, we find that doped impu- 
rities lead to a reduction of the spin gap in the dimer liquid 
phase e.g. for A = 0.5. In the VBC phase the effect is much 
less pronounced. We believe that the existence of an impurity- 
spinon bound state at A = 0.5 (see above) is responsible for 
this trend: if spinons get bound to available impurities (see 
below) the energy cost to create two of them is reduced. 



C. Spin density modulations 

We have also examined the spinon density map (equivalent 
to the local magnetization map) in the two spinon-two im- 
purity doped = 48 cluster, again for the same values A 
as above. The features shown in Figs. [TO|a-c) are consistent 
with our previous results for a single impurity. In the A = 0.5 
dimer liquid, a spinon clearly forms a bound state around each 
impurirty. In contrast, in Figs. 10 b) & (c) which we might 



view as some local region of a Valence Bond Glass , spinons 
seem to be repelled (at short distance) from the impurities. At 
further distances, strong spin density modulations are induced 
by a strong (average) dimerization (which reversely is weakly 
affected by the small magnetization of the sample). 



V. APPLICATION TO NMR EXPERIMENTS 

A. General considerations 

The different magnetization maps in the liquid and VBC 
phases are expected to lead to different characteristic NMR 
(theoretical) spectra which can be confronted to experiments 
in Herbertsmithite'^. The theoretical Copper NMR spectra 
is defined by the histogram of the local site magnetizations. 
Since Oxygen is located between two coppers, the effective 
local magnetization is obtained by summing up the contri- 



7 



butions from both sites (named here as bond magnetization). 
Both site and bond magnetization histograms have been com- 
puted for A = 0.5, A = 0.9 and A = 1 using the spinon 
density maps of Figs. [Tj^a-c) (one impurity on 48 sites) and 
Figs. 10 'a) & (c) (two impurities on 48 sites) providing the 
theoretical Copper and Oxygen NMR low-temperature spec- 
tra shown in Figs 1 1 Note that, for convenience, we work 
here at fixed magnetization'*" corresponding to one spinon per 
impurity i.e. m — ^riimp. This is a natural choice at very 
small temperatures if spinons bind to impurities or if the spin 
gap is very small. In any case, this corresponds physically to 
a magnetic field h > As{h — 0). For convenience, the Cop- 
per and Oxygen Knight shifts (x-axis) have been normalized 
w.r.t. the bulk values corresponding to a uniform spatial dis- 
tribution of the magnetization (mimicking high temperatures). 
Note however that, even for such a uniform distribution of the 
Copper spin density, the Oxygen nucleus located on the bonds 
connected to the impurity site see only one Copper atom in- 
stead of two, yielding a small "satellite" peak in the spectrum 
of Fig. [TTJb) at half the bulk Oxygen NMR shift. 

Comparing quantitatively high and low -temperature (in fact 
here T — 0) spectra obtained at a low fixed magnetic field 
would in principle require the knowledge of the tempera- 
ture dependent magnetic susceptibility x{T) which is beyond 
our numerical capability. However, a number of robust fea- 
tures of these theoretical NMR spectra can be linked to phys- 
ical behaviors of the spinons around impurities. E.g. when 
spinons spontaneously bind to impurities (as it is the case in 
the A = 0.5 dimer liquid phase within our model), one expects 
low-energy magnetic excitations, possibly at arbitrary small 
energy scale so that the broad theoretical spectra (at fixed m) 
of Fig. 1 1 c,d,i,j) should provide a good sketch of low-field 
NMR spectra. On the other hand. Figs. [TT|e,f) correspond 
to the opposite case where spinons are strongly delocalized 
around impurities leading to narrower structures in the spec- 
tra (in the same units). However, if the system retains a finite 
spin gap (as expected in this type of dimer liquid), for a small 
magnetic field, the average position of the spectrum should be 
strongly pushed downwards w.r.t. its high temperature refer- 
ence values. 



B. Herberthsmithite described as a Valence Bond Glass 

Let us now discuss our results in the light of the recently 
obtained ^^O NMR experimental spectra on the Kagome 
QHAF Herberthsmithite material which contains an intrinsic 
small Zinc substitution in the Copper Kagome planes'^. Al- 
though the experimental NMR spectra are strongly broaden 
by quadrupolar effects typical of powders, two separate struc- 
tures can clearly be followed from high to low temperatures. 
The average shift of the "satellite" peak is roughly half of the 
main one so that this peak has been interpreted as resulting 
from the resonance of the four Oxygen nucleus (labelled as 
"a" in Fig. 
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next to an impurity site. Our results show 
that dimerization is enhanced by impurity doping so that the 
proposed Valence Bond Glass^^ should be stabilized in an en- 
larged interval of A compared to the region of stability of the 
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FIG. 1 1 : (color online) Copper (left) and Oxygen (right) NMR spec- 
tra computed from the magnetization maps of Fig.|7|a-c) obtained for 
a single impurity on 48 sites (effective concentration of ~ 2%) and 
of Fig.[TO[a) and (c) obtained for two impurities on 48 sites (effective 
concentration of ~ 4%). The Knight shift variables are normalized 
w.r.t. the values corresponding to an ideal uniform distribution of the 
spin density and defining the reference NMR spectra (a,b) labeled 
as "high-temperature". Zero-temperature spectra for ~ 2% impurity 
concentration A = 0.5 (c,d), A = 0.9 (c,d) and A = 1 (g,h); Same 
for ~ 4% impurity concentration A = 0.5 (i,j) and A = 1 (k,l). 
For a single impurity, the impurity peaks are labelled according to 
the positions of the resonating nucleus shown in the insets (see also 
Fig.Hc)). 



VBC in the pure model. 

It has been recently argued that, confined spinons in the 
Valence Bond Glass can form a random-singlet phase with a 
wide distribution of spin gaps down to zero energy, ensuring 
the presence of fluctuating spins at arbitrarily low tempera- 
tures as indeed observed experimentally. The low-temperature 
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low-field magnetic response of a Valence Bond Glass is there- 
fore dominated by isolated impurities in rare regions of the 
sample, each of them carrying an almost free spinon loosely 
bound by a (shallow) VBC confining potential^'. Within this 
scenario the T, /i — > NMR spectrum is qualitatively given 
by the dominant magnetic response of the above regions. If 
this physics is at play in Herberthsmithite, Fig. [TT|h) should 
therefore give a crude qualitative sketch of the T ^ Oxy- 
gen NMR spectrum. This spectrum vaguely shows two broad 
structures (as experiments) but the assignment of the small 
shift structure to the a Oxygen is not so clear as shown in 
Fig.[TT|h). Note that the NMR spectrum would become more 
complicated when raising temperature since more and more 
regions of the sample will contribute (as e.g. Fig.[TO|c) giving 



Fig. 111)). 



VI. DISCUSSIONS AND CONCLUSIONS 

Theoretically, our approach is certainly based on an approx- 
imate framework relying on i) the truncation within the NN 
VB basis, ii) an overlap expansion up to a finite order and iii) a 
"minimal" spinon effective hamiltonian. However, we believe 
it still captures a number of important features and unambigu- 
ously shows that even a small (experimentally unavoidable) 
impurity doping should have important experimental conse- 



quences. Fig. 10 c) probably gives a fairly reliable picture of 



what is to be expected in a small area of the experimental sys- 
tem at finite magnetization. Even an impurity concentration 
as small as 4% leads to an important increase of the average 
dimerization and to a fairly inhomogeneous spin density dis- 
tribution. This property should be robust in a more elaborate 
treatment of the NN QHAF, like for example when using a 
refined spinon Hamiltonian including longer-range hoppings. 
It has been recently argued that, confined spinons in the 



VBC can form a random-singlet phase with a wide distribu- 
tion of spin gaps^^^, ensuring the presence of fluctuating spins 
at arbitrarily low temperatures as indeed observed experimen- 
tally. Interestingly, our data on small clusters show that the 
VBC pattern is not strongly confining, probably due to the 
weakness of the VBC amplitude compared to the spinon ki- 
netic energy. On the scale of the typical impurity separation, 
spinons tend to delocalize around impurities. However, dif- 
ferent types of dimerization patterns get pinned around dif- 
ferent impurities. In this sense, our findings have some simi- 
larity with the picture of the Valence Bond Glass^^*. Not sur- 
prisingly, we have not found any appreciable reduction of the 
spin gap (compared to the pure case) for the (two-) impurity 
configuration we have chosen (except when spinon-impurity 
bound states form at A = 0.5) : we believe our cluster is still 
too small to realize (rare) configurations with isolated impu- 
rities that could exhibit (almost) free spin within some (large) 
confining length. This issue certainly needs further work. 

On the experimental side, a more realistic description of 
Herbertsmithite may require to go beyond the simple NN 
QHAF. It is possible that longer range exchange interactions 
destabilize the VBC ground state in favor of the Z2 dimer liq- 
uid or some other candidate GS. It might also happen that the 
VBC phase is only stable at temperatures lower than those ex- 
perimentally reached. Finally, we note that, experimentally, a 
large enough DM anisotropjEl could fill out the spin gap of 
the Z2 dimer liquid or VBC phases. 
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